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Assuming a phenomenological self-energy 7mE(u) ~ \ui\ , ((3 — 1), which becomes gapped below 
T c , we derived a new gap equation. The new gap equation contains the effect of the kinetic energy 
gain upon developing a superconducting order parameter. However, this new kinetic energy gain 
mechanism works only for a repulsive pairing potential leading to a s-wave state. In this case, 
compared to the usual potential energy gain in the superconducting state as in the BCS gap equation, 
the kinetic energy gain is more effective to easily achieve a high critical temperature T c , since it is 
naturally Fermi energy scale. In view of the experimental evidences of the d-wave pairing state in 
the hole-doped copper-oxide high-T c superconductors, we discuss the implications of our results. 
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Numerous transport measurements of the high-Tc su- 
perconductors fjj indicate a very strong quasi-particle 
scattering rate (1/t(u>)) up to ~ 2000cm -1 or even higher 
up at the normal state, which is definitely electronic 
origin and is a source of most of the anomalous nor- 
mal state properties such as ac/dc conductivity, NMR, 
Knight shift, ARPES, Raman spectra, etc S. This 
unique quasi-particle scattering rate is a finger print of a 
poor metallic state due to the nearness to the Mott in- 
sulator. On the other hand, this strong scattering rate is 
gapped below the superconducting critical temperature 
T c or some higher temperature (T* > T c ) for most of 
the underdoped compounds, as indicated by all of the 
experiments above mentioned. Therefore, there is a clear 
indication that the quasi-particle scattering rate is sup- 
pressed as a kind of gap develops in the electronic channel 
and this gap is closely related to the superconducting gap 

On the other hand, many theories have been proposed 
to account for such high critical temperature in the high- 
Tc cuprates. Certainly, it is beyond the scope of this 
paper to judge the level of success of these theories. One 
of the common ingredients of many theories is that such 
a strong correlation effect, as indicated in the many nor- 
mal state properties, should somehow be involved in the 
superconducting pairing mechanism. Along this line of 
thinking, perhaps the Interlayer Pair Tunneling (IPT) 
theory, proposed early on by P.W. Anderson et. al. @, 
most actively exploited this idea to the limit, such that 
this theory requires a spin-charge separated non-Fermi 
liquid ground state (called a Luttinger liquid). As a re- 
sult, the c-axis normal transport is dynamically blocked 
(named as confinement) but only a momentum conserv- 
ing pair tunneling process is allowed via virtual process 
H . More intuitive account of this theory is that the frus- 
trated kinetic energy along the c-axis in the normal state, 
because of its Luttinger liquid ground state, is liberated 
in the superconducting state and this kinetic energy gain 
is the main source of the high critical temperature ||. 
In this theory, however, there are some of the issues yet 
to be cleared out, such as a realization of the Luttinger 
liquid ground state, for example. Nevertheless, the idea 
of the kinetic energy gain as a source of superconducting 
condensation energy seems to be supported with various 
experiments, in particular, such as the c-axis ac conduc- 
tivity and the c-axis penetration depth measurements [Q . 

In this paper, we take a phenomenological approach 
to the problem how to incorporate the kinetic energy 
gain upon entering a superconducting phase in the su- 
perconducting pairing mechanism, not necessarily con- 
fining in the c-axis kinetic energy. We assume that the 
kinetic energy frustration of the system is represented by 
a self-energy 7mE(w) ~ M' 3 , and this self-energy be- 
comes gapped upon developing a superconducting order 
parameter. We considered mainly (3 — 1 case since it is 
the most relevant case for high-T c superconductors ■ 



We obtained a new gap equation as a natural extension 
of the BCS gap equation when the condensation energy is 
drained not only from the potential energy but also from 
the kinetic energy as possibly realized with the above as- 
sumption. Oddly enough, the new gap equation shows 
that the kinetic energy gain works only for a repulsive 
pair potential and best for an isotropic s-wave state. For 
an attractive pair potential it works destructively and 
with a strong anisotropy of the order parameter the ef- 
fect quickly disappears. We show some numerical results 
of Tc as a function of the strength of kinetic energy frus- 
tration both for a s-wave and a d-wave state in a model 
Tc equation. The results indicate possible competition 
between a s-wave and a d-wave state in high-T c super- 
conductors. We discuss the implications of our results in 
view of experiments. 

We start with a usual BCS Hamiltonian. 
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where and c are the electron creation and annihila- 
tion operators, is the electron energy including chem- 
ical potential, V is a two body pair potential. Now, as- 
suming U(l) symmetry breaking order parameter, A = 
v 11k';\z , \<u D < c fc'aC-fc',T >i and introducing the 

Nambu spinor, %)A — (c£ j,c_y), we can readily obtain 
the free energy in the Matsubara frequency summation 
as follows. 
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where r = UbT and w„ = 7rr(2n + 1) is the Matsubara 
frequencies. 

Now, minimizing the above free energy with respect to 
A, i.e., SF/SA =0, we obtain the usual BCS gap equation 
as follows. 
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where E\ = ((% + A* A). 

This is a standard procedure to obtain the self- 
consistent mean field equation for the order parameter 
presumed for a given Hamiltonian. We want now to in- 
clude the frustrated kinetic energy and its regaining effect 
in the superconducting state in the above derivation of 
the gap equation. We do this by including the quasi par- 
ticle renormalization parameter Z(u ni A) in the above 
free energy and minimizing it with respect to A with 
special attention to the A dependence of Z(uj ni A) . We 
assume the form of a self-energy as follows, which is both 
experimentally and theoretically motivated ||. 
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ImYj(uo) — a\u>\ for 3A < \u\ < oj c 

for \lu\ < 3 A (4) 

where lu c is an ultraviolet cutoff. For a d-wave gap, more 
reasonable one would be like JmE(u) ~ oj 4 for \ui\ < 3A . 
However, even in that case, we can use the above simplifi- 
cation because it makes no qualitative change in our main 
results ||. The self-energy in the Matsubara frequency, 
E(w„, A), is obtained by the spectral representation, 



£(w„,A) = -- 
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and uj n is accordingly renormalized. 
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With this u)„ in the free energy, Eq.(2), we get 
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Finally, we obtain the gap equation for a s-wave order 
parameter, A = const. 
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The Tc equation is obtained by taking A — > limit. 
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where ifi = Z n (A — 0). This is our main result and 
now let's analyze the above gap equation. The term, 
^18Zn\ in the numerator in the right-hand side of Eq. 
(9) is the new term originating from the kinetic energy 
gain; otherwise, Eq. (9) is a simple BCS Tc equation 
with just an ordinary self-energy correction, Z n , in the 
denominator. As well known, if not the new term, the 
above equation can have a solution only for an attrac- 
tive interaction (i.e. V < 0). However, with the new 
term, it can have a solution even with a purely repulsive 
interaction (i.e. V > 0) since — 18 ~ 2.9 with a ~ 0.5 
(which is a reasonable value for the high-T c cuprates ) 
and Z n > 1 always. What if V < then ? It is a dis- 
aster. The new kinetic energy gain term works only for 
a repulsive pair potential but it works destructively for 
an attractive pair potential. Therefore I think that the 
newly found self-energy derived term contains some more 
physics than the kinetic energy gain, as discussed later. 



In the case of d-wave order parameter, Hi nt and A{k) 
is defined as follows. 
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Here, V(k, k') is a positive definite pair potential in mo- 
mentum space and A(fc) is a sign-changing order param- 
eter of d-wave symmetry. The gap equation can be de- 
rived similarly by taking a variation of the free energy, 
Fo-wave, with respect to A(fc); it is slightly more in- 
volved than the s-wave case since the order parameter 
has a structure. The result is written as, 



A(fc) 



EE 



n k 



V{k,k')A{k') 



V(k,k )7 M+(3Arii „ )a ] 5A 

-t 2^ 2^ — — ^ — A 



n k' 



+ El,] 



SA(k' 



(11) 



Note that the second term contains now f^9p) , which 
has just measure zero since it equals <5fe^ ax ,fc' sgn{A{k')); 
it would be just 1 in the s-wave case. Therefore, the 
new kinetic energy gain term practically has no contri- 
bution in the d-wave gap equation; more generally, its 
effect quickly vanishes for any anisotropic gap A(k) un- 
less A max has a finite support in k— space [[To). Now in 
order to gain some more insight about the newly found 
term, we recall the self-consistency definition of the or- 
der parameter, A(fc) =< ^ k ,V{k,k')ck> [C-k>] >■ At 
first sight, the second term in Eq. (11) appears to be 
a redundant term. Although it is not rigorously proved 
here, this extra term doesn't ruin the self consistency but 
actually completes it. We note that the right-hand side 
of Eq. (11) is the pair susceptibility ( a kind of response 
function ) with respect to the U(l) symmetry breaking 
source, A(fc). The first term in the right-hand side of 
Eq. (11) is the usual pair susceptibility including the 
self-energy correction, Z(oj n ,A). It is well known that 
a vertex correction is necessary whenever a self-energy 
correction is included in order to satisfy some continu- 
ity equation (or underlying symmetry ) for any response 
function. Therefore, we can consider the self-energy de- 
rived second term as a corresponding vertex correction 
- it indeed looks like so. And in this sense, the new 
extra term should contain not only the kinetic energy 
gain effect but also some other complicate effect (general 
backflow). 

Now let us show some numerical results of a model 
Tc equation. All calculations are performed in two di- 
mensional momentum space assuming a circular Fermi 
surface for simplicity of calculations. In Fig.(l) we show 
Tc as a function of a, the strength of inelastic scattering, 
both for s-wave and d-wave states with different values 
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of the coupling constant, A, (A = VN(0), N(0) is the 
density of states at Fermi level). For all calculations, the 
ultraviolet cutoff scale uj c is taken to be O.beV and the 
BCS potential cutoff u>d to be 0.5eV since the source 
of the pairing interactions must be an electronic origin. 
This choice of parameters is mainly for exemplary pur- 
pose, though. For a s-wave state, the results of A = 0.1 
and 0.2 are shown - extremely weak repulsive potentials, 
otherwise we would get much too high T c . The message is 
clear; we can easily achieve several hundreds K of T c with 
a very weak repulsive interaction and it increases with 
more kinetic energy frustration (i.e., with a larger value 
of a) as expected from the Tc equation, Eq. (9). For a 
comparison, we show a d-wave case (solid square) . For a 
d-wave state, we drop the second term in the gap equa- 
tion, Eq. (11), as explained above and the Tc equation is 
obtained by taking A(k) — ► limit. For simplicity of nu- 
merical calculations, we assume V(k, k') — V\sin(<f) — <f)')\ 
and A(k) = A max cos(2((>) in two dimensional momen- 
tum space with a circular Fermi surface, and (f> is an 
angle along the Fermi surface. With uj c = O.beV and 
u>d = 0-beV as in the s-wave case, we need A = 1.5 
(an order of magnitude stronger interaction) in order to 
get a comparable Tc (~ 100i<Q for a ~ 0.4 - 0.6. This 
value of a is the experimental value near optimal dop- 
ing with Tc ~ 90K in YBCO @. Tc quickly decreases 
with a increasing since now a just enters the ordinary 
pair-breaking self-energy correction in contrast to the s- 
wave case. The results show that if the ground state of 
the system is a marginal Fermi liquid (or even when only 
some part of Fermi surface is marginal Fermi liquid like) 
there is a possible competition between a s-wave state 
and a d-wave state for the experimentally relevant range 
of a, provided there is a weak repulsive interaction in the 
s-wave channel [|ll) . 

Some remarks and speculations are in order. First, in 
this paper we considered only (3=1 case (7mS(oj) ~ 
\o-il 13 )- we can easily extend our analysis to the other 
power. In general, SZ n (A) /SA ~ A 13 . As a result, if the 
system is a Fermi liquid ((3 = 2), the kinetic energy gain 
term simply drops out of the Tc equation |)J. On the 
other hand, if (3 is any power less than 1 , the kinetic en- 
ergy gain term becomes singular as A — > and the gap 
equation becomes a non-linear equation. In this sense, 
(3=1 case is marginal not only for breaking Fermi liquid 
Pj but also for the superconducting instability. Second, 
it is straightforward to extend our formalism to include 
c-axis dynamics by introducing an anisotropy in and 
a cylindrical shape of order parameter A(k). The whole 
formalism goes as the same as in two dimension. There is 
no kinetic energy gain effect for a d-wave order parameter 
in regard to enhancing Tc as the system undergoes from 
the incoherent state to the coherent state. The inter- 
layer pair tunneling mechanism by Anderson et. al. [|]|^] 
is certainly a different approach from ours. It is not clear 
how to reconcile two approaches at present. Third, as 



a pure speculation and viewing our numerical results, it 
may be that the s-wave instability and the d-wave insta- 
bility closely compete in the marginal Fermi liquid state 
depending on the potential strength of the s-wave and 
d-wave channels. Then our results might be relevant to 
the electron-doped high-T c compound |Q. 

In conclusion, assuming a phenomenological self- 
energy, 7mE(w) ~ |oj |^ , ((3 = 1), we derived a new gap 
equation which includes the kinetic energy gain effect. It 
is argued that the new gap equation is a natural extension 
of the BCS gap equation with a consistent vertex correc- 
tion when the self-energy is a functional of the gap func- 
tion. For the corresponding Tc equation, it is also shown 
that the new kinetic energy gain term becomes relevant 
only when (3 < 1, and (3 = 1 is the marginal case. When 
(3=1, the new Tc equation can provide surprisingly high 
critical temperature for the case of an isotropic s-wave or- 
der parameter with a pure repulsive potential. For the 
d-wave case, the new kinetic-energy gain term drops from 
the Tc equation because of the anisotropy of the order 
parameter. Numerical results indicate that there can be 
a competition between the s-wave and d-wave instabil- 
ities in a marginal Fermi liquid state depending on the 
strength of pair potential in s-wave and d-wave channels. 
In view of the facts that there is a substantial strength 
of repulsive interaction in s-wave channel besides the in- 
teraction in d-wave channel in the Hubbard or t-J model 
and its relevance to high-T c superconductors, our results 
may have some relevance to the electron-doped high-T c 
compound. Considering the unexpected behaviors of the 
new gap equation, it is desirable to have further investi- 
gations of it. 
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is irrelevant for the Tc equation as far as 7 > 1 in 
ImT,(uj) ~ \lu\~ 1 for \u>\ < 3A since the effect of any higher 
power term (7 > 1) drops out in the Tc equation as tak- 
ing A — + limit. However, it makes a difference for the 
determination of gap value A(T) from the gap equation. 
3^(1°) m Eq.(ll) is correct only for the Tc equation 
because of the same reason as explained in [9]. In gen- 
eral, we need a functional derivative, —fj^jpy^-- Then 
A(fe) < A max would contribute to the detailed shape 
of JmE(oj) for \u>\ < 3A ma x and quite generally we can 
guess JmE(aj) ~ |w| 7 with 7 > 1 for any anisotropic or- 
der parameter. Then as argued in [9], only A max part 
matters for the Tc equation. A possible exception is 
the case when A max spans a finite region in momentum 
space. 

Some might argue that the s-wave channel should be at- 
tractive in real metals if the phonon mediated interaction 
and the retardation are taken into account. The retarda- 
tion is not compulsory, though. If the system can lower 
the energy by developing superconducting order param- 
eter without taking advantage of the retardation, it will 
do so. 

Most of experiments on the electron-doped compound 
(NCCO) indicate a s-wave state; however, the normal 
state properties of it show more Fermi liquid like behav- 
iors (private communication with Jae H. Kim). 
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FIG. 1. Tc vs a = a\u\). For all data, 

ljd = 0.5eV and u c = 0.5eV; solid square (d-wave, A = 1.5), 
open circle (s-wave, A = 0.2), and open diamond (s-wave, 
A = 0.1). 
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